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We present a possible extension of the usual relativistic nuclear mean field models widely used to describe 
nuclear matter towards accounting for the influence of possible intrinsic fluctuations caused by the environ- 
ment. Rather than individually identifying their particular causes we concentrate on the fact that such eff'ects 
can be summarily incorporated in the changing of the statistical background used, from the usual (extensive) 
Boltzman-Gibbs one to the nonextensive taken in the form proposed by Tsallis with a dimensionless nonexten- 
sivity parameter q responsible for the above mentioned effects (for g — > 1 one recovers the usual BG case). We 
illustrate this proposition on the example of the QCD-based Nambu - Jona-Lasinio (NJL) model of a many-body 
field theory describing the behavior of strongly interacting matter presenting its nonextensive version. We check 
the sensitivity of the usual NJL model to a departure from the BG scenario expressed by the value of \q - 1|, in 
particular in the vicinity of critical points. 

PACS numbers: 



I. INTRODUCTION 



BG distribution, f(E), are connected as follows: 



In all studies of relativistic propeities of nuclear matter, 
mean field models are usually the models of first choice [[lllltl- 
These models use a statistical approach based on Boltzmann- 
Gibbs (BG) statistics which is, strictly speaking, only correct 
when the corresponding heat bath is homogeneous and infi- 
nite. These conditions are by no means met in realistic situ- 
ations in which nuclear matter occurs. Usually one encoun- 
ters some inherent problems arising, for example, from the 
smallness of the collisional systems and their rapid evolution. 
These, among other things, render the spatial configuration of 
the system far from uniform and prevent global equilibrium 
from being established (cf. |3] and references therein). As 
a result, some quantities become non extensive and develop 
power-law tailed rather than exponential distributions. The 
widely used way to account for these effects is to resort to 
a nonextensive statistics, known as <7-statistics [4]. The new 
phenomenological nonextensivity parameter q occuiiing there 
is supposed to account for all possible dynamical factors vio- 
lating the assumptions of the usual BG statistics. This is re- 
covered in the limit of ^ — > 1 . Because it enters into the re- 
spective formulas of the particular dynamical model used for 
a given investigation, it allows for a simple phenomenological 
check of the stability of the model against possible deviations 
from the BG approach. 

So far, applications of the nonextensive approach are nu- 
merous and cover all branches of physics Q. These in- 
clude high energy multiparticle production processes (cf., |0]) 
and different aspects of nuclear and quark matter (l7l-|9l. The 
nonextensive framework can also be derived from a special 
treatment of kinetic theory investigating complex systems in 
their nonequilibrium stationary states [ lOJ.^ome examples of 
more specialized topics can be found in jl ill and references 
therein. For an illustration, the Tsallis distribution, hq{E), and 



exp, 



/(£) = iexp(-|) 



1-(1 



(1) 

(2) 



It is usually argued that, for the ^ > 1 case, q - \ \s a. mea- 
sure of intrinsic fluctuations of the temperature in the system 
considered LlZl . whereas ^ < 1 is usually attributed to some 
specific correlations limiting the available phase space ifisll or 
to the possible fractality of the allowed phase space iflill (other 
possible inteipretations were considered in [11])'. 

For our further considerations of importance are recent ap- 
plications of nonextensive statistics in description of nuclear 
01811 and quarkonic matter [8, 19, 20], the later of which we 
shall continue here. In lilSll . the ^-version of the Walecka 
many-body field theory flj] has been investigated. It was 
shown there that <7-statistics results in the enhancement of the 
scalar and vector meson fields in nuclear matter, in dimin- 
ishing of the nucleon effective mass and in hardening of the 
nuclear equation of state (only the ^ > 1 case was consid- 
ered there). In [8^ the relativistic equation of state of hadronic 
matter and a quark-gluon plasma at finite temperature and 
baryon density was investigated in the framework of nonex- 
tensive statistical mechanics. In our work [19] we investi- 
gated a nonextensive version of another mean field theory, 
namely the QCD-based Nambu - Jona-Lasinio (NJL) model of 
a many-body field theory describing the behavior of strongly 
interacting matter presented recently in ll2l]| . This time, unlike 
in [|l8[l . we used the quark rather than the hadronic degrees 
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' One must admit at this point that this approach is subjected to a rather hot 
debate of whether it is consistent with the equilibrium thermodynamics or 
else it is only a handy way to a phenomenological description of some in- 
trinsic fluctuations in the system fl^ . It should be therefore noticed that it 
was demonstrated on general grounds iUtIi that fluctuation phenomena can 
be incorporated into a traditional presentation of a thermodynamic. The 
Tsallis distribution l[T] belongs to the class of general admissible distribu- 
tions which satisfy thermodynamical consistency conditions and which are 
therefore a natural extension of the usual BG canonical distribution l|2). 
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of freedom and, because of this, we had to consider both the 
q > I and q < I cases. This g'-NJL model allowed us to dis- 
cuss the ^-dependence of the chiral phase transition in dense 
quark matter, in particular the quark condensates and the ef- 
fective quark masses and their influence on the masses of n 
and cr mesons and on the spinodal decomposition (cf., lfl9ll 
for details). These results helped us proceed further and con- 
sider critical phenomena in strongly interaction matter using 
^-statistics (these phenomena are of interest nowadays, cf., 
for example, ll23il24 ll, but were so far not investigated in non- 
equilibrium environment provided by <7-statistics). In partic- 
ular, we shall now concentrate on the influence of dynamical 
factors causing nonextensivity and represented by the param- 
eter q in the vicinity of the critical end point (CEP). 



II. BASIC ELEMENTS OF THE g-NJL MODEL 



First we present the basic elements of the g'-NJL model in- 
troduced in ifigl l (to which we refer for more details). 



A. The usual NJL model 



We start with the usual QCD based NJL model based on 
BG statistics discussed in |21]. It is the standard SU{3) NJL 
model with U{1)a symmetry described in, with the usual La- 
grangian of the NJL model used in a form suitable for the 
bosonization procedure (with four quarks interactions only), 
from which we obtain the gap equations for the constituent 
quark masses M,: 



Ml 



(3) 



with cyclic permutation of j, k - u,d,s and with the quark 
condensates given by (qiqi) = -iTT[S ,(/?)] {S i(p) is the quark 
Green function); m, denotes the current mass of quark of fla- 
vor /. We consider a system of volume V, temperature T 
and the quark chemical potential yu, characterized by the 
baryonic thermodynamic potential of the grand canonical en- 
semble (with quark density equal to p, = Ni/V, the baryonic 
chemical potential /^b = ^(/z„ +i^d+ fJ-s) and the baryonic mat- 
ter density as pg = \(pu+Pd+ Ps)), 



^l(T,V,fii) ^E-TS 



(4) 



The internal energy, E, the entropy, S , and the particle num- 
ber, Ni, are given by ^ (here = ^Mf + p^): 

E = — 2^ I p dp (l-n, -n,) - 

-gsV ^ {{qiqi)f -2gDV{uu){dd){ss), (5) 

i-ii,d,s 

S = ~V 2 f p^dp-S, (6) 

where S = [n,Tnn; -H (1 - n,)ln(l - «,)] + 
+ [ni 1 - n,]. 



Ni 



Nc r 



p dp (ni - n,) . 



(7) 



The quark and antiquark occupation numbers ai^e, respec- 
tively. 
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and with them one calculates values of the quark condensates 
present in Eq. (|3]l. 



i=uM,s 



dp. (10) 



Eqs. (O and (fTOl ) form a self consistent set of equations from 
which one gets the effective quark masses M, and values of 
the corresponding quark condensates. 

The values of the pressure, P, and the energy density, e, are 
defined as: 



with 



n(//,-, r) 



E(jXi,T) 
£(lJi,T)^ — - — (11) 
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P(0,0) = e(0,0) = 0. 



B. The q extension of the NJL model - the ^-NJL 

The <7-statistics is introduced by using the ^-form of quan- 
tum distributions for fermions (+1) and bosons (-1) in Eqs. 
(O and (|9]l. This is done following a prescription provided in 
III5L namely by replacing n and h by 
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(12) 



" ~eg(/3{Egi-tii))±V 
(the important point to notice is that one encounters here E^i = 

sj^gi P^^ ^^'^^ because of M^i also energy is now a q- 
dependent quantity). Denoting x - /3(E - fi) one has that for 

q > I 



eJx) 



{\ + {q - \)x\i-^ if x>0 



[\+{\ -q)x]~ if x<Q 
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whereas for q < I 



III. RESULTS 



[I + (q - l)x]i-' if jc < 
e,(x) = <j _ . (14) 

[l+(l-^)x]T^ if x>Q 

This is because only then can one consistently treat on the 
same footing quarks and antiquarks (and for all values of x). 
This should show the particle-hole symmetry observed in the 
^-Fermi distribution in plasma containing both particles and 
antiparticles, namely that 



ng(E,/3,iJ.,q) = 1 - n2-q(-E,/3, -ju). 



(15) 



This means, therefore, that in a system containing both parti- 
cles and antiparticles (as in our case) both q and 2 - q occur 
(or, when expressed by a single q only, one can encounter both 
^ > 1 and ^ < 1 at the same time). These dual possibilities 
warn us that not only q > I but also q < 1 (or (2-q) > 1 have 
physical meaning in the systems we are considering. This dif- 
ferentiates our£-NJL model from the ^-version of the model 
presented in fW\. Notice that for ^ — > 1 one recovers the 
standard FD distribution, n(fi, T). Actually, it is important to 
realize that for T — > one always gets «^(yu, T) — > n{fi, T), 
irrespectively of the value of ^ [18]. This means that we can 
expect any nonextensive signature only for high enough tem- 
peratures (how high depends on circumstances and on the kind 
of observable considered, for illustration of this point see re- 
sults presented in our paper ifiol and Fig. |2]below). 

Our ^-NJL model is then obtained by replacing the formu- 
las of Section UlAl with their ^-counterparts in what concerns 
the form of the FD distributions. Additionally, when calcu- 
lating energies and condensates we follow [7, 8] and use the 
^'-versions quark condensates, replacing Eqs. dTHt . ([3) and (|5]l 
by their q'-forms: 



i=u,a,s 



I 



(1 



dp, (16) 



Mqi = mi - 2g^{qiqi)^ - 2g ^{q jq j) ^{q^qt) ^ 
+ niiM, 
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If,..'. 



-9sV 2 {(^i9i\f-2gDV(uu\{dd)^(-ss\^. (18) 



-u,d,s 



On the other hand, again following Qlltl, densities which are 
given by the the g'-version of Eq. (|7]l are calculated with riq's 
(not with rig, as in (fTST l and in (fTSIl). The pressure for given q 
is calculated using the above Eg and the ^-entropy version of 
Eq. dSl) with (cf. ill) 



^1 ^ 1"? + (1 - ^qif ln?(l - "<?'■)] + 



l-rin 



(19) 



Eq. ( |T6] l together with the ^-version of the gap equation, Eq. 
([TTT i, are the basic equations from which one deduces all re- 
sults presented here. 



Before presenting our results concerning nonextensive criti- 
cal effects we shortly repeat the previous results (cf., lll9ll20l ). 
In Fig. [T]we present the typical pressure at critical tempera- 
ture Tcr obtained in a ^-NJL model as a function of compres- 
sion pIpq calculated for different values of the nonextensivity 
parameter q (see [19] for more details on spinodial decompo- 
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FIG. 1 : The pressure at critical temperature r„ as a function of com- 
pression pi Pa calculated for different values of the nonextensivity 
parameter q (the area marked at the upper panel is shown in detail in 
the lower panel). The dots indicate positions of the inflection points 
for which first derivative of pressure by compression vanishes. As in 
[21] for ^ = 1 the corresponding compression is p/po = 1.67 (and 
this leads top = 318.5 MeV); it remains the same fox q> 1 consid- 
ered here (but now p = 321 MeV for ^ = 1.01 and p = 326.1 MeV 
for q = 1.02) whereas it is shifted to p/po = 1-72 for q < I (ji = 313 
MeV for q = 0.99 and p = 307.7 MeV for q = 0.98). 

sition and chiral symmetry restoration in ^-NJL model)^. No- 
tice that the effect is stronger for q < I and that, essentially. 



^ There is still an ongoing discussion on the meaning of the temperature in 
nonextensive systems. However, in our case the small values of the param- 



4 



the saddle point remains at the same value of compression. 
When one moves away from the critical temperature, the typ- 
ical spinodal structure occurs, which is more pronounced for 
lower temperatures whereas its sensitivity to the q parameter 
gets stronger with increasing temperature (cf., [19il ). How- 
ever, it occurs that, for each temperature (even for very small 
one) a ^ > 1 exists for which there is no more mixed phase 
and for which the spinodal effect vanishes. This seems to be a 
quite natural effect in the scenario in which ^ > 1 is attributed 
to the fluctuations of the temperature in a system considered 
as proposed in [12]. On the contrary, effects like correlations 
or limitations of the phase space considered in lil3i ll4] work 
towards an increase of the Tcr and make the spinodal effect 
more pronounced. 

A few remarks are in order here (for more detailed discus- 
sion we refer to lfl9ll ). Nonextensive dynamics enter the NJL 
calculations through the quark (antiquark) number distribu- 
tion functions Uq, ifiqi). These functions are connected with 
the respective quark (antiquarks) spectral functions in the NJL 
model. However, deviations from the exponential shape of q- 
exponents, as defined in Eqs. ( fT3] l and (fT4l i. are negligible for 
values of q close to unity (in our case 0.98 < q < 1.02). It is 
also important to notice that Eqs. (fT3T l and (fT4l i are symmetric 
for q \ -q. The differences between q < \ and q> \ cases 
observed in our results are then due to our way of defining the 
energy ( fTSl l and entropy ([19), which, following [L,8J, we do 
by using rf'qi and h'^'qi instead of iiqi and hqj ^ . Because now 
for q < \ distributions if'qj and h'^'qi are closer to unity than 
Hqj and hqi, therefore the absolute values of quark condensates 
(as given by Eq. (fTSl l) begin to decrease for q - 0.98 at lower 
temperature as compared with the q - \ case. Therefore the 
corresponding energy is larger, which means that ^ < 1 in- 
troduces some residual attractive correlations which rise the 
energy and lead to hadronization occurring at lower tempera- 
ture. On the other hand, q > \ introduces fluctuations which 
decrease the effective occupations {n''qi and h^'qi) and the en- 
ergy, and smears out the chiral phase transition. In Fig. |2] 
we present our phase diagram in the -T plane for different 
nonextensivity parameters considered here with positions of 
the corresponding critical end points (CEP) for different val- 
ues of q clearly indicated. The overlap of curves observed in 
Fig. |2] (inlet) indicates how the critical end point is smeared 
to a kind of critical area. This is because fireballs created in 
different events can have different values of q (representing, as 
mentioned before, action of all factors responsible for the de- 
parture of our system from the usual BG approach - not spec- 
ified here in detail but, in general, resulting in specific corre- 



eter q deduced from the data allow us to argue that, to first approximation, 
= r used here and in 1 18]. In high energy physics it is just the hadroniz- 
ing temperature (and instead of the state of equilibrium one deals there with 
some kind of stationary state). For a thorough discussion of the temperature 
of nonextensive systems, see | 22]. 
^ It is worth notice that in |18], which considers only the q > I case and 
uses number distributions without powers of q, the significant elfects were 
obtained only for much larger values of the nonextensive parameter q = 
1.2. 



lations of quarks or fluctuations of temperature mentioned be- 
fore). Therefore when analyzing experimental data one most 
probably will encounter such a critical area instead of a well 
defined CEP. 
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FIG. 2: Phase diagram in the i^-NJL model in T - fi plane for values 
of q considered before: q = 0.98, 1.0 , 1.02. Solid and dashed lines 
denote, respectively, first order and crossover phase transitions. The 
results are presented for three different values of the nonextensivity 
parameter q with the vicinity of the ((^-dependent) critical end points 
(CEP) enlarged in the inlet. The crossover phase transition for q = 
0.98 and for yu — > takes place for a smaller temperature T. 

The role of all these factors is shown in more detail in Fig. |3] 
which shows the baryon compression p/po (calculated in the 
vicinity of the critical values of temperature and density in- 
dicated by the corresponding dotted lines) as the function of 
the chemical potential p for different values of the nonexten- 
sivity parameter, q - 0.98, 1.00, 1.02. Notice the remarkable 
difference of the density derivative at the critical point: from 
the smooth transition through the critical point for ^ < 1 to 
a big jump in density for critical value of chemical potential 
for q > I. It reflects the infinite values of the baryon number 
susceptibility, ^g: 



Xb 



i=u,d,s 



2 (£1-2 




(20) 



The transition between confined and deconfined phases and/or 
chiral phase transition ||23| can be seen by measuring, event 
by event, the difference in the magnitude of local fluctuation 
of the net baryon number in a heavy ion collision [25]. They 
are initiated and driven mainly by the quark number fluctua- 
tion, described here by xb, and can survive through the free- 
zout [25 1. Consequently, ourq-NJL model allows us to make 
the fine tuning for the magnitude of baryon number fluctua- 
tions (measured, for example, by the charge fluctuations of 
protons) and to find the value of the parameter q characteristic 
for this system. However, it does not allow us to differentiate 
between possible dynamical mechanisms of baryon fluctua- 
tion. We close by noticing that using q dependent leads 
to ^-dependent parameter e of the critical exponents which 
describe the behavior of baryon number susceptibilities near 
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FIG. 3: The baryon compression p/po (calculated in the vicinity of the critical values of temperature and density indicated by the coiTesponding 
dotted lines) as function of the chemical potential ji for different values of the nonextensivity parameter, q = 0.98, 1.00, 1.02. The summary 
presented in the top-left panel is detailed in the three consecutive panels. 



the critical point ll26ll . Whereas in the mean field universal- 
ity class one has e = e' = 2/3, our preliminary results us- 
ing the g'-NJL model show a smaller value of this parame- 
ter for ^ > 1, (e ~ 0.6 for q=1.02) and greater for q < I 
(e ~ 0.8 for q=0.98). It would be interesting to deduce the 
corresponding values of q from different models and compare 
them with results on a lattice which, by definition, should cor- 
respond to ^ = 1 (it should be mentioned at this point that 
there are already attempts to apply Monte Carlo methods, sim- 
ulating lattice gauge field dynamics as based on non-extensive 
rather than extensive thermodynamics, see li27ll and references 
therein). 

In order to further investigate the q dependence ofxB let us 
rewrite Eq. ( |20l l in the following form (recall that p = Ng/V 
and A^^ is <7-version of Eq. (Q): 

The q'-versions of occupation numbers, n^, and n^,-, are taken 
from Eq. ( fT2] i. The ^-version of energies there depend on 



masses M^,-, which are given by gap equation ( [TtI i in a quite 
involved way. Therefore, the ^-dependence enters here in two 
ways: by rather straightforward replacement of exp(...) by 
the respective exp(...) in Eq. ( fT2l i and by quite involved q- 
dependence of M^, given by the gap equation ([TtI i. Therefore, 

1 r - n 

Xb(Mb,T) = -T^f \X (Mb) + X (Mb)\ (22) 

with 
where 

f(q) = (2-q) if (q - 1)(E,, - i^b) > 0, 

fiq) - q. otherwise. 
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FIG. 4: Upper panel: the chemical potential (/jg) dependence of the 
light quarks condensate in the vicinity of the critical region calcu- 
lated according to Eq. ( I16t for different values of the nonextensivity 
parameter q: q = 0.99, 1.0 and 1.01. Bottom panel: the jig depen- 
dence of the chemical potential derivative of the light quark mass 
M,„ calculated according to Eq. Ml\ in the critical region for the 
same values of q as above. 



Our results are presented in Figs. 2] and |5] It turns out that 
the chiral phase transition investigated here (Fig. |5) is mainly 
driven by the behavior of the light quark mass derivative, see 
Fig. m which in turn is determined by the behavior of the 
light condensate, cf.. Fig. [3] Thus the dynamic of the nonex- 
tensive effects is generated not so much by the nonextensive 
form of occupation numbers in Eq. (fTSl i but rather by the main 
gap equation (O where both the condensates and the effective 
quark masses are present. 



rv. SUMMARY 



FIG. 5: The dependence of the baryon number susceptibility, Xb, 
in the vicinity of the critical region, calculated according to Eq. i2\\ 
for different nonextensivity nonextensive parameters, q = 0.99, 1.0 
and 1.01. Notice that it is essentially identical with results presented 
in the bottom panel of Fig. |4] 



required by the application of the BG approach by using the 
TsalUs version of nonextensive statistical mechanics [4]. All 
factors causing this departure are summarily described by the 
nonextensivity parameter q, such that q-l quantifies departure 
from the BG situation (which is recovered for ^ ^ 1). 

We have investigated two possible scenarios correspond- 
ing to ^ > 1 and q < I, respectively, which, as mentioned, 
coiTespond to different physical interpretations of the nonex- 
tensivity parameter For q < 1 (usually connected with 
some specific correlations |ll3|] or with fractal character of the 
phase space 1 14]) we observe a decreasing of pressure, which 
reaches negative values for a broad (<7-dependent) range of 
temperatures and increasing of the critical temperature The 
^ > 1 case (usually connected with some specific nonstatis- 
tical fluctuations existing in the system |12J) we observe a 
decreasing of the critical temperature, Tcrit, and therefore in 
the limit of large q we do not have a mixed phase but rather 
a quark gas in the deconfined phase above the critical line 
(on the contrary, the compression at critical temperature does 
not depend on q. As in 1 18] the resulting equation of state is 
stiffer (in the sense that for a given density we get larger pres- 
sure with increasing q). As expected, the effects depend on 
the temperature, and tend to vanish when the temperature ap- 
proaches zero. Fig. |3]shows that the nonequilibrium statistics 
dilutes the border between the crossover and the first order 
transition. Finally, Figs. |4] and |5] demonstrate that the most 
important <7-dependence is coming from the main gap equa- 
tion (O, where both the condensates and the effective quark 
masses are present, rather than from the the nonextensive form 



We have investigated the sensitivity of critical behavior of 
the QCD based NJL type of mean theory type presented in 
ll2lll . the g'-NJL model, to the departure from the conditions 



It acts therefore in the same way as including of the Polyakov loop into the 
NJL model [2M. 
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of occupation numbers in Eq. (fT2] l. 

We would like to end by stressing that our results could 
be of interest for investigations aimed at finding the critical 
point in high energy heavy ion collisions [3] or when studying 
part icularities of the equation of state (EoS) of compact stars 
029(1 ■ The fact that they depend on the parameter q means that 
the exact position of such a point or the type or shape of EoS 
could be quite different from what is naively expected. 
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